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ABSTRACT 

The objective of this work was to investigate the damage mechanisms in composite bonded 
skin/stringer constructions under uniaxial and biaxial (in-plane/out-of-plane) loading conditions as 
typically experienced by aircraft crown fuselage panels. The specimens for all tests were identical 
and consisted of a tapered composite flange, representing a stringer or frame, bonded onto a 
composite skin. Tests were performed under monotonic loading conditions in tension, three-point 
bending, and combined tension/bending to evaluate the debonding mechanisms between the skin 
and the bonded stringer. For combined tension/bending testing, a unique servohydraulic load 
frame was used that was capable of applying both in-plane tension and out-of-plane bending loads 
simultaneously. Specimen edges were examined on the microscope to document the damage 
occurrence and to identify typical damage patterns. The observations showed that, for all three load 
cases, failure initiated in the flange, near the flange tip, causing the flange to almost fully debond 
from the skin. 

A two-dimensional plane-strain finite element model was developed to analyze the different 
test cases using a geometrically nonlinear solution. For all three loading conditions, principal 
stresses exceeded the transverse strength of the material in the flange area. Additionally, 
delaminations of various lengths were simulated in two locations where delaminations were 
observed. The analyses showed that unstable delamination propagation is likely to occur in one 
location at the loads corresponding to matrix ply crack initiation for all three load cases. However, 
the current two-dimensional plane-sham finite element model may not fully account for the 



complex three-dimensional damage pattern observed. A detailed investigation of this damage 
pattern may require a local three-dimensional analysis of the damaged area. 
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INTRODUCTION 

Carbon epoxy composite structures are widely used by today's aircraft manufacturers to 
reduce weight. Co-curing, co-bonding, and secondary bonding have become the most promising 
processes to replace traditional mechanical fastening methods. Composite materials have been 
introduced fairly recently into primary structures of commercial airplanes. Tire failure processes in 
composites are not as well understood as in metals. Previous investigations of the failure of 
secondary bonded structures were conducted with specimens cnt from a full-size panel to verify the 
integrity of the bondline between the skin and the flange or frame [1 ] . However, these panels were 
rather expensive to produce and there was a need for a test configuration that would allow detailed 
observations of the failure mechanism at the skin/flange interface. A simpler specimen 
configuration was proposed in reference 2 with the investigations focusing on the failure 
mechanisms of a bonded skin/flange configuration loaded in bending and transverse shear [2-5]. 
These loading conditions may be appropriate for a variety of applications, but in many cases 
composite structures may experience both bending and membrane loads during in-flight service. 
The failure mechanisms under multi-axial loading may be complex in that they do not represent a 
simple combination of the various load components but involve interaction between the loads. 

The first objective of this work was to investigate the damage mechanisms in composite 
bonded skin/stringer structures under monotonic tension, three-point bending, and combined 
tension/bending loading conditions. For combined tension/bending testing, a unique 
servohydraulic load frame was used that was capable of applying axial tension and transverse 
bending loads simultaneously [6, 7]. Microscopic investigations of the specimen edges were 
performed to document the damage occurrence and to identify typical damage patterns. 

The second objective of this work was to develop an analytical methodology to accurately 
predict the damage observed during the experiments. All three load cases were analyzed using a 
detailed two-dimensional plane-strain finite element model. Both linear - and geometrically nonlinear 
simulations were performed. A stress analysis was used to predict the location and orientation of 



the first transverse crack based on the principal transverse tension stress distribution in the flange 
tip area. A fracture mechanics approach was utilized to determine when a delamination would grow 
from this transverse crack. Mode I and mode II strain energy release rate contributions, G : and 
G n , were calculated for all load cases using the virtual crack closure technique [8, 9] and compared 
to existing mixed-mode fracture toughness data. 


EXPERIMENTAL INVESTIGATION 
MATERIALS AND SPECIMEN PREPARATION 

The specimens tested in this investigation consisted of a bonded skin and flange assembly 
as shown in Figure 1. All specimens were machined from the same panels and were similar to the 
specimens used in the previous monotonic and fatigue tests reported in references 2 through 5. 
Both the skin and the flange laminates had a multidirectional lay-up. The skin lay-up, consisting of 
14 plies, was [0/45/90/-45/45/-45/0] s and the flange lay-up, consisting of 10 plies, was [45/90/- 
45/0/90] s . 

Both skin and flange were made from IM6/3501-6 graphite/epoxy prepreg tape with a 
nominal ply thickness of 0.188 mm. First, the flange and skin laminates were cured separately. 
The flange parts were then cut into 50.0 mm long strips. A previous investigation had indicated 
that the angle at the flange tip has a significant influence on the strength [2], A 20° taper angle was 
suggested in a previous investigation, however, a 27° taper was machined along the edges. 
Subsequently, the flange was adhesively bonded to the skin using a 1 77 °C cure film adhesive 
from American Cyanamid (CYTEC 1515). A grade-5 film was used to yield a nominally 
0.1 27 - mm thick bondline. However, because some of the adhesive flowed outwards during cure, 
the bondline thickness averaged 0.102 min. A diamond saw was used to cut the panels into 
25.4-mm wide by 203.2-mm long specimens (the specimen length is the only difference to 
specimens used in previous studies [2,3,5] which were 127.0 mm long). Each specimen was then 
equipped with two strain gages, one located in the center of the flange and the other located on the 
skin as close to the flange tip as possible (Figure 1). Ply properties and adhesive material 
properties were measured at Boeing and are part of the standard design database for the V-22 tilt- 
rotor aircraft. Typical material properties for the composite tape and the adhesive material as taken 
from reference 2 are summarized hi Table 1. 



EXPERIMENTAL PROCEDURE 


A total of five tension tests were performed in a servohydraulic load frame in displacement 
control. The actuator speed was controlled at 0.4 mm/min. The specimens were mounted in 
hydraulic grips to give a gage length of 127.0 mm as shown in Figure 2. An extensometer with a 
25.4 mm gage length was mounted on the backside of the specimen and centered on the flange tip 
as shown in Figure 2. The tests were terminated when the flange debonded unstably from one of 
the flange tips. 

Five bending tests were performed in a servohydraulic load frame at a monotonic rate of 
1.52 mm/min. A photograph of the three-point bending test fixture is shown in Figure 3. The 
configuration used was similar to the one used in previous studies [2, 5]. The bottom support had 
a 127.0 mm span. Mid-span deflections were recorded using a spring loaded direct current 
displacement transducer (DCDT) contacting the center of the flange as shown in Figure 3. The 
tests were stopped after the flange debonded unstably to the center of the specimen. 

Tests under combined membrane and bending loading conditions were performed in the 
axial tension and bending (ATB) servohydraulic load frame shown in Figure 4. In this ATB load 
frame designed at the NASA Langley Research Center, the axial load cell is incorporated in the top 
grip that rotates with the upper specimen part (see close-up in Figure 5) [6, 7], The specimens 
were mounted into the machine with great care placed on correct alignment of the specimen and the 
top grip/load cell set-up. The specimens were initially preloaded in load control to an axial tension 
load of 17.8 kN, which is 85% of the average damage initiation load determined for the tension 
test. While maintaining this preload, a transverse bending load was then applied in displacement 
control until flange debonding occurred. Maximum specimen deflections at the top grip contact 
point were recorded using a spring loaded linear variable differential transformer (LVDT). 

The first specimen was tested with a gage length of 127.0 mm. The transverse load was 
applied at a constant rate of 1.52 mm/min (in accordance with three-point bending tests). As the 
stroke was increased, the specimen failed near the lower grip without evidence of flange 
debonding. As a result, the gage length was reduced to 101.6mm for the four remaining 
specimens. The shortening allowed the lower flange tip to be closer to the bottom grip which 
increased the bending moment at the lower flange tip so that flange debonding occurred before skin 
failure. The transverse displacement rate was tripled to 4.57 mm/min to reduce testing time caused 
by the large specimen deflections necessary for damage initiation. 



TEST RESULTS 


In Figures 6 to 11, typical results of each test are shown as plots of load versus 
displacement or strain versus load. For each loading configuration, the plots of load versus 
displacement and strain are from the same specimen. The strain-load curves are shown for flange 
strain and skin strain. The loads, flange and skin strains are reported in Tables 2 and 3 (tension 
and three-point bending tests) for the point of possible damage initiation as well as for the flange 
debond load. In Table 4 (ATB tests), the results are presented for the flange debond load only, 
since no damage initiation prior to failure was observed in these tests. 

For tension specimens, the load-displacement curves were slightly nonlinear over a wide 
range as shown in Figure 6. Possible damage initiation was assumed when a small initial load 
drop was observed prior to flange debonding. At this point, a crack in one flange tip or even a 
small delamination along one flange comer was observed, hi one specimen, no initial load drop or 
visible damage could be detected prior to flange debond. In general, the initial load drop occurred 
above 90% of the flange debond load. Figure 7 shows the slightly nonlinear strain-load response 
until flange debonding. In all specimens, a load drop was also accompanied by a decrease in strain. 

Typical plots for the three-point bending tests are shown in Figures 8 and 9. The load 
versus mid- span deflection curves of all three- point bending specimens showed nonlinear behavior 
at higher loads indicating possible damage initiation. Sometimes this behavior was accompanied by 
a minor load decrease (see Figure 8). However, no cracks or delaminations could be observed 
prior to ultimate flange failure. In all specimens, the nonlinearity or initial load decrease was again 
detected above 90% of the flange debond load. Both skin and flange strains showed linear 
behavior before flange debonding (see Figure 9). No decrease in load or strain was observed on 
the strain-load response prior to flange debonding. 

In contrast to tension and three-point bending tests, the transverse load versus transverse 
displacement curves obtained from the ATB tests showed no indication of damage formation until 
just prior to skin failure for all specimens. Flange debonding could not be identified from these 
plots. A typical example is shown in Figure 10. Figure 11 shows a characteristic strain versus 
transverse load response starting from the strain caused by the axial load. Although Figure 11 
shows linear behavior, the strain curves in some specimens deviated slightly from linearity prior to 
flange debonding. Flange debonding as detected during the experiments always occurred at 
maximum flange strain and was sometimes accompanied by a drop in skin strain. 



MICROSCOPIC INVESTIGATION 


Photographs of the polished specimen edges were taken under a light microscope after 
testing was completed to document the occurrence of matrix cracks and delaminations. Damage 
was documented based on location at each of the four flange corners identified in Figure 12,_ 
Comers 1 and 4 and corners 2 and 3 had identical damage patterns. Typical damage patterns which 
were similar for all three loading configurations are shown in Figure 13. These drawings are 
based on the micrographs taken after the tests. Figure 14 presents two such photomicrographs for 
a three-point bending specimen. In general, failure in tension and three-point bending specimens 
occurred on one flange tip only, with no clear preference for either flange tip. Due to the moment 
distribution resulting from the loading and boundary conditions of the ATB test, failure in the ATB 
specimens occurred at the flange tip with the higher bending moment only, i.e., the flange tip 
closer to the lower grip. 

At comers 1 and 4, a delamination miming in the 90745° flange ply interface (delamination 
A) initiated from a matrix crack in the 90° flange ply as shown in Figure 13(a). At longer 
delamination lengths, new matrix cracks formed and branched into both the lower 45° as well as 
the upper 90° flange ply. However, no branching into the adhesive bondline was observed. 

At corners 2 and 3 a matrix crack formed at the flange tip hi the 90° flange ply that 
subsequently ran through the lower 45° flange ply and the bondline into the skin. Subsequently, a 
split (delamination Bl) formed from the tip of that matrix crack within the top 0° skin ply as 
depicted in Figure 13(b). In some cases, a second delamination (delamination B2) was observed 
below the first in the top 0°/45° skin ply interface. Both delaminations were present over a long 
distance until delamination Bl stopped and delamination B2 continued. 


ANALYTICAL INVESTIGATION 
FINITE ELEMENT ANALYSIS 

The finite element (FE) method was used to analyze the test specimens for each loading 
case. The goal of this investigation was to study damage initiation using a stress analysis and the 
potential for delamination propagation using a fracture mechanics approach. FE models for one 
undamaged and two damaged specimens were developed and loads and boundary conditions were 
applied to simulate the three load cases. The two-dimensional cross section of the specimens was 
modeled using eight-noded quadrilateral plane strain elements using quadratic shape functions and 
a reduced (2x2) integration scheme. The models are shown in Figures 15 and 16. For the entire 
investigation, the ABAQUS geometric nonlinear analysis procedure was used. For the tension and 



three-point bending loading cases, the results of linear analyses were compared to those of 
nonlinear analyses. For the ATB test, only the ABAQUS nonlinear solution was used since this 
allowed the axial load to rotate with the specimen as it deformed under the transverse load and 
accounted for the membrane stiffening effect caused by the axial load. 

For the model of the undamaged specimen, a refined mesh was used in the critical area of 
the 90° flange ply where cracking was observed in the test specimens. An outline and two detailed 
views of the mesh are shown in Figure 15. Outside the mesh refinement area, all plies were 
modeled with one element through the ply thickness. In the refined region, two elements were used 
per ply thickness except for the first three individual flange plies above the bondline and the skin 
ply below the bondline which were modeled with four elements. Three elements through the 
thickness were used for the adhesive film. The model consisted of 6492 elements and 19975 nodes 
and had 39888 degrees of freedom. The properties used to simulate the behavior of the 
graphite/epoxy material and the adhesive were measured at Boeing and are part of the standard 
design database for the V-22 tilt-rotor aircraft. Typical properties as taken from reference 2 are 
summarized in Table 1 . 

Based upon the experimental observations shown in Figures 13 and 14, a "damaged" 
model was also developed that included discrete matrix cracks and delaminations. The mesh 
described for the undamaged specimen was also used for this model, except for the critical area 
around the flange tip where delaminations were modeled as shown in Figure 16 (a) for corners 1 
and 4 and Figure 16 (b) for corners 2 and 3. The initial matrix crack was modeled perpendicular to 
the flange taper. Damage was modeled at one flange tip only as shown in Figure 16. At the 
opposite taper, the mesh used in the model of the undamaged specimen was employed. This 
procedure was used to simulate the occurrence of damage onset only. In order to keep tire FE mesh 
simple and avoid skewed elements, the split in the 0° ply (marked as "Delamination B 1 " in Figures 
13 (b) and 14 (b)) was modeled as a delamination propagating between the adhesive film and the 
top 0° ply of the skin (Figure 16 (b)). It is assumed that this slight alteration in geometry does not 
significantly influence the computed energy release rates. For the prediction of delamination 
growth, the fracture toughness values obtained from standard specimens with unidirectional layup 
were used as described below. It is inherent to a two dimensional plane strain FE model that the 
geometry, boundary conditions and other properties are constant through the entire width. This 
may not always capture the true nature of the problem. As shown in Figures 13 and 14, the 
delamination pattern changed from comer 3 to comer 4 from a delamination running in the 90745° 
interface to a delamination propagating between the adhesive film and the top 0° ply of the skin. 
This three dimensional effect can not be accounted for in the current model. 

The schematics of the specimen, boundary conditions, and load cases applied in the 
simulations are shown in Figure 17 for the tension and three-point bending cases and in Figure 18 



for the combined tension and bending case. For the tension and bending case, the mean loads 
reported for the point of damage initiation in Tables 2 and 3 were applied. At this point, matrix 
cracks are likely to form. For the simulation of the combined tension and bending loads in the ATB 
test, the top grip, the load cell, and the load pin were also modeled using three-noded quadratic 
beam elements as shown in Figure 18. A rectangular beam cross section was selected to model the 
square cross section of the top grip and load pin and a circular beam cross section was used for the 
model of the cylindrical load cell. The beams were connected to the two-dimensional plane strain 
model of the specimen using multi-point constraints to enforce appropriate translations and 
rotations. As shown in Figure 18, nodes 1-29 along the edge of the plane strain model were 
constrained to move as a plane with the same rotation as beam node A. To be consistent with the 
ATB tests, a constant axial load, P, was applied in a first load step while transverse loads remained 
zero, hi a second load step, the axial load was kept constant while the load orientation rotated with 
the specimen as it deformed under the transverse load. In the FE simulation this transverse load 
was applied as a prescribed displacement which corresponded to the mean of the transverse stroke 
(31 mm) reported in Table 4. For the beam model of the steel parts (top grip, load cell, and load 
pin), a Young's Modulus of 210 GPa and a Poisson's Ratio of 0.3 were used as material input 
data. 

The Virtual Crack Closure Technique (VCCT) described in references 8 and 9 was used to 
calculate strain energy release rates for the delaminations. The mode I and mode II components of 
the strain energy release rate, G t and G n , were calculated as 


and 
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where A a is the length of the elements at the crack tip, X/ and Y) are the forces at the delamination 
tip at node i, and u m ’ and v m ' are the relative displacements at the corresponding node m behind the 
delamination tip as shown in Figure 19. Similar definitions are applicable for the forces at node j 
and displacements at node 1. Both forces and displacements were transformed into a local 
coordinate system (x', y'), that defined the normal and tangential coordinate directions at the 
delamination tip in the deformed configuration. The total strain energy release rate, G r , was 
obtained by Slimming the individual mode components as 

G t = Gj +G n . 


( 3 ) 



The mode III component is zero for the plane strain case. The data required to perform the 
VCCT in equations (1) to (3) were accessed from the ABAQUS result file. The calculations were 
performed in a separate post processing step using nodal displacements and nodal forces at the 
local elements in the vicinity of the delamination front. 


ANALYSIS RESULTS 
Global Response 

First, the global response of the specimens was computed at load levels corresponding to 
the damage initiation loads observed in the experiments. The load-displacement and the load-strain 
behavior computed for all three load cases were compared to the corresponding experimental 
results. This global response was used to examine whether the FE model, the boundary 
conditions, the load cases and the material properties used in the model were accurate. Note that the 
experimental data only represent one typical specimen, thus not accounting for any experimental 
scatter. Displacements were reported at the locations where they were taken in the experiments. 
Strains were computed at a single location corresponding to the center of the strain gage. A 
schematic of the deformed geometries, the boundary conditions, and the load cases applied in the 
simulations are shown in Figure 20 for all three load cases. 

In the schematic of the deformed FE tension model in Figure 20(a), the elongation of the 
specimen caused by the applied tensile load is shown along with the bending effect caused by the 
load eccentricity. The load versus displacement plot in Figure 21 shows that the linear and 
nonlinear FE simulations are in good agreement. Moreover, there is very marginal difference 
between the analyses and the experiments, the model therefore accurately captured the global 
response. In Figure 22, a comparison of measured strains and computed results is shown. The 
strain-load responses for the skin are again in good agreement between the experiments and both 
analyses. For the flange strains the nonlinear analysis yielded results which were in agreement with 
the experiments. Strains obtained from the linear analysis were veiy low. Hence, a geometrically 
nonlinear FE analysis is necessary to account for the effects due to the load eccentricity in the 
flange region, the asymmetric layup with respect to the neutral axis and the membrane stiffening 
effect. 

In the three-point bending test, the vertical displacement at the flange tip (Figure 3) was of 
the order of the skin thickness for the load level investigated. Hence, a geometrically nonlinear FE 
solution procedure may also be needed in this case. Both linear and nonlinear analyses were 
performed and computed displacements and strains were compared. Tire load-displacement plot in 
Figure 23 and the strain-load plot in Figure 24 show that both simulations and the experiments are 



within 10% of each other. Since the linear and nonlinear analyses are within 1%, a linear analysis 
is sufficient for the load level investigated, i.e., the mean load obtained form experiments for 
damage initiation. Higher load levels, however, might require a nonlinear analysis. 

For the tension/bending loading case the load-displacement and strain-load plots in Figures 
25 and 26 show that the nonlinear FE simulation and the experiments are within 20% of each 
other. A possible explanation for the stiffer response of the numerical model based on the 
difference in axial modulus measured in tension and bending is given in the appendix. 


Local Response 

Tire local response was studied in the critical area of the 90° flange ply where cracking was 
observed in the experiments. The goal of these FE analyses was to investigate damage initiation 
using a stress analysis and the potential for delamination propagation using a fracture mechanics 
approach. 

The stress analysis was used to study the initial failure in the form of matrix cracks from 
which delaminations may start to grow. As shown in Figures 13 and 14, the first crack always 
occurred in the 90° flange ply closest to the skin. Results from previous investigations [2, 3] 
suggested that the local force and moment resultants, may be used to characterize damage onset. As 
shown in Figure 27, force and moment resultants per unit width were calculated from the normal 
stress a xx at the flange tip. Resultants were computed for damage initiation and flange debond load 
level based on stresses obtained from linear and nonlinear FE analyses. Results plotted in Figure 
28 show that the three-point bending test is free of axial tension as the force resultants /V xx are zero 
as expected. A bending moment is present in the tension specimen caused by the load eccentricity 
in the flange region and the asymmetric layup of the combined skin and flange laminate with 
respect to the neutral axis. For the ATB test, calculated force and moment resultants lie between 
computed pure tension and pure bending values. 

Earlier investigations [2, 3] also indicated that the maximum ply principal transverse tensile 
stress may cause the initial failure hi the form of matrix cracks from which delaminations may stall 
to grow. Failure, therefore, may occur when the computed principal tensile stress in the 2-3 plane 
of the ply (normal to the fiber direction) exceeds the transverse tensile strength of this ply. In this 
model, the ply 2-3 plane for a 90° ply coincides with the global x-y plane of the model. Maximum 
principal stresses can therefore be taken straight from the finite element results. The vector plot in 
Figure 29 shows the trajectories of the maximum principle tensile stresses in the flange ply. At the 
-45790° and 90745° interface, multiple vectors are displayed since the stresses were not averaged 
across boundaries for elements with different material properties. Comparing the trajectories in 



Figure 29 with the damage patterns shown in Figures 13 and 14 shows that the crack starts to 
grow perpendicular to the trajectory of the maximum principle tension stress. Computed maximum 
principal tension stresses a max in the elements with labeled element numbers in the 90° ply in 
Figure 29 are presented in Table 5 and compared in Figure 30. For all three loading conditions, 
maximum principal tensile stresses computed for the damage initiation load level have similar 
magnitudes. Towards the center of the ply, principal stresses exceeded the failure strength of 
61.1 MPa as found for a similar type of material (AS4/3501-6 in [10]) and hence ply cracks may 
develop. Consequently, the stress analysis based on the comparison of computed stresses with 
failure strengths appeal's to be an appropriate method to determine the location of the initial failure 
and the orientation of the resulting crack. 

A fracture mechanics approach was used to investigate delamination growth once the initial 
crack had formed. It had been shown that care must be exercised in interpreting the values for G : 
and G n obtained using the virtual crack closure technique for interfacial cracks between two 
orthotropic solids [11]. Therefore, it had been suggested to use element lengths A a at the crack tip 
in such a manner that the computed results are insensitive to the variation of A a For practical 
applications the element size (length and height) should not be less than 1/10 of a ply thickness 
because the assumption that each ply can be modeled as a continuum is no longer valid. The 
element size at the crack tip also should not exceed the ply thickness as this requires smearing 
properties of individual plies. For the current investigation, the element length A <7 was chosen to be 
1/4 of the ply thickness for the delamination in the 90745° flange ply interface and 1/3 of the 
bondline thickness for the simulated propagation along the 0° skin ply/adhesive fihn interface. 

The initial crack was modeled on one flange tip perpendicular to the flange taper as 
suggested by the microscopic investigation and the stress analysis described in the previous 
paragraph. Recall that the models of the discrete cracks and delaminations are shown in Figure 16. 
During the analyses, the delaminations were extended by adding new nodes at the crack tip and in 
front of the crack tip. These nodes were then assigned to the elements on one side of the crack thus 
creating a row of disconnected elements which simulated the delamination. Strain energy release 
rates were computed at each tip location for the flange debond load observed in the experiment. 
The delamination lengths, a, are measured from the end of the initial matrix crack as shown in 
Figure 16. For comers 1 and 4 (delamination in the 90745° flange ply interface), the delamination 
was extended to 0.6 mm which corresponds to a length where matrix crack branches were 
observed in the experiments as shown in Figures 13(a) and 14(a). The results are plotted in 
Figures 31 through 33 for all three loading conditions. Initial mode I and mode II values, 
computed at delamination onset (a = 0.034 mm, i.e., delamination length equal to the length of 
the first element as shown in Figure 16(a)), are similar for each type of test as shown in Table 6. 



In Figures 31 through 33, G n for all load cases increases monotonically while G, begins to level 
off. Computed mixed mode ratios G[/G T for all the load cases are shown in Figure 34. 

For the simulated propagation along the 0° skin ply/adhesive film interface (comers 2 and 
3) the delamination was extended to 1.6 mm. Plots of strain energy release rates computed for 
various delamination lengths are given in Figures 35 through 37. The computed values at 
delamination onset (a= 0.04 mm, i.e., delamination length equal to the length of the first element 
as shown in Figure 16(b)) are given in Table 7. Comparing all load cases, computed G l results 
appear to have similar magnitudes. Mode II values, however, differ noticeably. As shown in 
Figures 35 through 37, at a < 0.2 mm the mode I contribution decreases with increasing 
delamination length. The mode II contribution and the total strain energy release rate, however, 
increase with increasing delamination length. For a > 0.2 mm, the mode I contribution and the 
total sham energy release rate increase slowly with increasing delamination length for the bending 
load case mode and decrease slowly for the pure tension and the combined tension/bending load 
case. For all three load cases, the mode II contribution continues to increase at a much slower rate 
with increasing delamination length. Corresponding mixed mode ratios G,/G T are shown in Figure 
38. The results discussed above will be used in a mixed-mode failure investigation to determine 
whether delamination onset and unstable propagation are possible at the applied loads where 
damage was observed in the experiments. 


MIXED-MODE FAILURE INVESTIGATION 


Accurate mixed-mode failure criteria are necessary for the prediction of delamination 
growth. A bilinear mixed-mode failure criterion was suggested in reference 12 for AS4/3501-6, a 
material similar to IM6/3501-6. The mixed-mode failure response was presented by plotting the 
mode I component of the mixed-mode fracture toughness versus the respective mode II 
component. A different approach to present the data was suggested in reference 13 where 
mixed-mode fracture toughness values, G c , were plotted versus the mixed mode ratio G : /G T (see 
Figure 39). When this ratio is zero, G c is simply the mode I fracture toughness, G Ic . Alternatively, 
G c becomes the mode II fracture toughness, G llc , when the mixed-mode ratio equals unity. An 
equation resulting from a least square regression cubic curve fit to these data defines the 
mixed-mode delamination fracture criterion for each ratio as: 


G = 75.3 + 214.7- 


— | - 70.5- 

v gJ 


( . 


-2- + 327.4- 

G. 


v G x 


(4) 



Hence, for a given mixed-mode ratio, growth is possible when the total mixed-mode strain energy 
release rate exceeds the critical value. 

In the current study, computed total strain energy release rates, G T , were compared to the 
critical value, G c , for the appropriate mixed mode ratio (G,/G T ) for each load case in order to 
determine the potential for delamination growth. Values calculated for delamination onset 
( a= 0.034 mm, i.e., delamination length equal to the length of the first element as shown in 
Figure 16(a) )_in the 90745° flange ply interface are below the fracture toughness data as shown in 
Table 6 and Figure 40. Consequently, onset is unlikely to occur at the load corresponding to 
damage initiation for all three tests. Propagation in the 0° skin ply/adhesive film interface, on the 
other hand, will occur in all three cases as the computed results for delamination onset 
(i a= 0.04 mm, i.e., delamination length equal to the length of the first element as shown in 
Figure 16(b))_are much higher than the reported fracture toughness values as shown in Table 7 and 
Figure 41. It is important to recall that in the experiments the failure was observed as a split in the 
0° ply as shown in Figures 13 (b) and 14 (b). This failure was modeled as a delamination 
propagating between the adhesive film and the top 0° ply of the skin as shown in Figure 16 (b) 
only in order to avoid skewed elements. Hence, for predicting growth of the split in the 0° ply it is 
appropriate to compare energy release rates computed for the 0° skin ply/adhesive film interface 
with fracture toughness values obtained from standard tests. Unstable propagation is likely since 
the calculated G T -values remain above the mixed-mode fracture toughnesses over the entire 
simulated length as shown in Figures 35 to 37. This assumption is confirmed by the experimental 
results of this study. 

The above findings suggest that once a matrix crack has formed, a delamination 
(delamination B1 from Figure 14) will also form and grow in an unstable manner between the 
adhesive film and the top 0° skin ply. The second delamination observed in the 90745° flange ply 
interface (delamination A from Figure 14) requires more energy to initiate than available at the load 
levels used in this FE analysis, i.e., the loads corresponding to possible matrix ply crack initiation. 
The energy for this second delamination may come from an increase in load or may be caused by 
an increase in G x due to the presence of the first delamination in the 0° skin ply/adhesive film 
interface. These two possibilities are studied in detail in the following paragraphs. The response of 
the numerical model with respect to the difference in axial moduli measured in tension and bending 
is discussed in the appendix at the end of the paper. This discussion includes the effect of a 
lowered axial modulus on computed strain energy release rates for specimens subjected to bending 
and combined tension/bending loading. 

First, strain energy release rates were calculated at the flange debond load observed in the 
experiment. Results for the second delamination in the 90745° flange ply interface are shown in 



Figure 42 and are included in Table 6. Strain energy release rates are higher, however, they still 
remain below the fracture toughness. This suggests that delamination growth is unlikely. 

Second, the possibility of an increase in G T due to the presence of the first delamination 
was investigated. On one flange side a delamination of 10 nun was simulated between the adhesive 
film and the top 0° skin ply as an initial disbond. This length was kept constant during the 
remainder of the investigation. At the opposite taper, the delamination in the 90745° flange ply 
interface was modeled as described in the previous paragraph. The model is shown in Figure 43. 
Again, the second delamination was extended to 0.6 mm which corresponds to a length where 
matrix crack branches were observed in the experiments. The total strain energy release rates 
computed at loads corresponding to possible damage initiation are plotted in Figure 44 for all three 
loading cases. The values computed previously for the condition without the assumed initial 
disbond (G T from Figures 31 through 33) were included in Figure 44 for comparison. For both 
conditions, G T increases monotonically with growing delamination length. For the tension and 
combined tension/bending load case computed values are smaller for the condition where 
delamination formation was assumed between the adhesive film and the top 0° skin ply. The 
computed values at delamination onset (a= 0.034 mm, i.e., delamination length equal to the 
length of the first element as seen in Figure 43) are shown in Figure 45 for delamination onset and 
flange debond load. Values for both load levels are again below the fracture toughness data, which 
suggests that delamination growth hi the 90745° flange ply interface is also unlikely in the presence 
of a 10 mm long delamination in the 0° skin ply/adhesive film interface. Based upon this study it 
can be concluded that delamination growth in the 90745° flange ply interface cannot be explained 
using strain energy release rates based on two dimensional plane strain FE analyses. The observed 
three-dimensional delamination pattern seen in the specimens could not be accounted for with the 
current model. A detailed and more accurate investigation should include a local three-dimensional 
analysis of the delaminated area. 


CONCLUDING REMARKS 

The damage mechanisms in composite bonded skin/stringer constructions under uniaxial 
and biaxial (in-plane/out-of-plane) loading conditions have been investigated using experimental 
and numerical approaches. Tests were performed under monotonic loading conditions in tension, 
three-point bending, and combined tension/bending to evaluate the debonding mechanisms 
between the skin and the bonded stringer or frame. Microscopic investigations of the specimen 
edges showed that all tests yielded similar damage patterns. For all three load cases, failure initiated 
in the flange, near the flange tip, causing part of the flange to fully debond from the skin. 



Based upon the experimental findings, a two-dimensional nonlinear plane-strain finite 
element (FE) analysis was performed using the ABAQUS FE code. For tension and three-point 
bending tests, linear and geometrically nonlinear simulations were performed. Overall, both 
computed results were found to be in good agreement with the experimental data. To account for 
the large displacements observed in the ATB tests, only the geometrically nonlinear analysis was 
performed. Again, the results were in good agreement with the experiments. 

A stress analysis was used to investigate the onset of failure. This approach showed that 
the location and orientation of the initial transverse ply crack in the flange are dependent on the 
stress distribution in the critical area near - the flange tip. For all three loading conditions, computed 
maximum principal tensile stresses were almost identical and exceeded the transverse tension 
strength of the material. A fracture mechanics approach was used to determine the potential for 
delamination growth from the initial transverse crack. In this approach, delaminations of various 
lengths originating from the transverse crack as observed in the experiments were simulated. 
Mode I and mode II strain energy release rate contributions were calculated for all load cases 
using the virtual crack closure technique. Computed total strain energy release rates were compared 
to critical values obtained from an existing mixed-mode failure criterion. Tire results suggest that 
once a matrix ply crack has initiated in the flange, a delamination will form and grow in an unstable 
manner between the adhesive film and the top 0° skin ply as observed in the micrographs. A 
second delamination located in the 90745° flange ply interface requires more energy to initiate than 
was computed to be available at loads corresponding to possible damage initiation. Strain energy 
release rates calculated at the flange debond load observed in the experiment were higher than the 
initiation values but still remained below the fracture toughness. Computed values also remained 
below the fracture toughness data in the presence of a 10 mm delamination modeled in the 0° skin 
ply/adhesive film interface. Consequently, delamination growth in the 90745° flange ply interface 
cannot be explained using strain energy release rates computed from a two-dimensional plane strain 
FE analysis. A detailed and more accurate investigation of the observed complex delamination 
pattern needs to include a local three-dimensional analysis of the damaged area. 


APPENDIX 

In the global response analysis, the load-displacement and the strain-load behavior 
computed for all three load cases were compared to the corresponding experimental results. For the 
tension specimen the load versus displacement plot (Figure 21) and the strain versus load plot 
(Figure 22) show that there is little difference between the nonlinear analyses and the experiment. 
This also holds for the strain-load responses of the skin and the flange as shown in Figure 22. For 



the three-point bending specimen, the load-displacement plot in Figure 23 and the strain-load plot 
in Figure 24 show that the simulations and the experiments are within 10% of each other. The 
plots in Figures 25 and 26 show that the nonlinear FE simulation and the ATB experiment are 
within 20% of each other. 

The slightly stiffer response of the numerical model may be explained by the fact that the 
material data used in the FE simulation originate from the literature. For a consistent simulation, 
material data should be taken from the batch of material that was used to manufacture the 
specimens. Furthermore, the figures only represent one specimen of each load case and do not 
include any experimental scatter. The stiffness difference between the model and the experiments, 
however, is more pronounced for the bending and ATB test than for tension loading. This could be 
caused by the fact that the axial modulus, E n , was measured in tension and not in bending. As 
reported in [14], axial moduli depend on the type of test method used. Moduli obtained from 
bending tests may be up to 20% lower than moduli from tensile tests. This suggests that tensile 
moduli should be used for the simulation of tensile dominated problems while flexural moduli 
should be used when modeling bending problems. Therefore, the effect of a 20% lower axial 
modulus on the computed global response and on the strain energy release rates of the specimens 
subjected to bending and combined tension/bending loads was studied in detail. 


EFFECT OF LOWERED AXIAL MODULUS ON THE COMPUTED GLOBAL BEHAVIOR 

The effect of a 20% lower axial modulus on the global response for the bending and 
combined tension/bending load case was studied first. The new load-displacement and the 
strain-load responses were computed and compared to the corresponding experimental and 
previously computed results. For the three-point bending test, the load versus displacement plot in 
Figure 46 shows that the nonlinear FE simulations for the two different moduli differ by about 
20% as expected. The results suggest that using the lower modulus yields better agreement with 
the experiments. In Figure 47, a comparison of measured strains and computed results is shown. 
The strain-load responses for both skin and flange differ about 20%. Results obtained from the 
analyses with the lower axial modulus yield a 10% higher compliance. Computed results taken 
from the previous analysis, however, are about 10% stiffer than the experiments. Based on a 
comparison of the strain results a decision as to which modulus should be used in a simulation is 
not possible. For the ATB test, the load-displacement plot in Figure 48 shows that the nonlinear 
FE simulations for both moduli are almost identical. The strain-load responses for the skin and 
flange are shown in Figure 49. At the end of the tensile preloading, computed strains are about 
20% higher than the values computed previously and the values measured in the tests. The slope 



during the bending phase of the test is almost identical for both computations and the experiment. 
Ideally, for the simulation of an ATB test, moduli obtained from tensile and flexural tests should be 
used for input. The tensile modulus would be used to compute the membrane stiffness terms of the 
element stiffness matrix while the flexural modulus would be used to determine the flexural 
stiffness terms. 


EFFECT OF LOWERED AXIAL MODULUS ON THE COMPUTED STRAIN ENERGY 
RELEASE RATES 

For a fracture mechanics investigation it is also essential to know how a 20% lower axial 
modulus will affect the computed strain energy release rates of the specimens subjected to bending 
and combined tension/bending loading. As in the earlier investigations, the delamination in the 
90745° flange ply interface was extended to 0.6 mm. For the simulated propagation along the 
0° skin ply/adhesive film interface the delamination was extended to 1 .6 mm. For comparison, 
computed strain energy release rates from the earlier studies were included in the plots. The results 
for the three-point bending loading case are plotted in Figures 50 and 51. As expected, mode I 
and mode II values were about 20% higher in those cases where a 20% lower axial modulus was 
used in the simulations. For the ATB load case, plots of computed strain energy release rates with 
increasing delamination length are given in Figures 52 and 53. Computed mode I and mode II 
values were typically up to 10% higher when a 20% lower modulus was used in the simulations. 
The increase in the mode I component of the strain energy release rate is less pronounced than in 
the mode II component. This non-uniform behavior may be explained by the fact that only the 
axial modulus, E n , was modified in the simulations. The computed total strain energy release rates 
at delamination onset (a = 0.034 mm, i.e., delamination length equal to the length of the first 
element) are shown in Figures 54 and 55. Values remain below the fracture toughness data, 
suggesting that delamination growth in the 90745° flange ply interface is unlikely even if the axial 
modulus was 20% lower. Propagation in the 0° skin ply/adhesive film interface will occur as the 
computed results are higher than the values computed previously and exceed the reported fracture 
toughness values. 

The study shows that computed strain energy release rates are about 20% higher if a 20% 
lower axial modulus is used in the simulation of structures predominantly subjected to bending. 
For a combined load case, the increase is less pronounced and computed strain energy release rates 
are about 10% higher. These results may be used to estimate the effect of a lowered axial modulus 
on calculated strain energy release rates without performing additional EE analyses. 
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TABLE 1. MATERIAL PROPERTIES. 


IM6/3501-6 Unidirectional Graphite/Epoxy Tape [3] 

E n = 144.7 GPa 

E 2 2 = 9.65 GPa 

£ 33 = 9.65 GPa 

v 12 = 0.30 

V 13 = 0.30 

V 23 — 0.45 

G 12 = 5.2 GPa 

Go = 5.2 GPa 

G 23 =3.4 GPa 


CYTEC 1515 Adhesive 


E = 1.72 GPa 

v =0.30 

(assumed isotropic) 


TABLE 2. RESULTS FOR TENSION TESTS. 


Specimen 

Damage 

initiation 

load 

kN 

Damage 
initiation 
flange strain 
pe 

Damage 
initiation 
skin strain 

|X£ 

Flange 

debond 

load 

kN 

Flange strain 
at flange 
debond load 
pe 

Skin strain 
at flange 
debond load 
pe 

2 

20.5 

1225 

5619 

22.8 

1348 

6084 

4 

21.8 

1419 

6312 

23.6 

1231 

6685 

6 

19.9 

1185 

5834 

23.1 

1187 

6599 

8 

20.9 

1300 

6051 

23.0 

1115 

6463 

10 

21.1 

1360 

6092 

21.1 

1360 

6092 

Mean 

20.9 

1298 

5982 

22.7 

1248 

6385 

Std. Dev. 

0.7 

96 

264 

0.9 

105 

282 

CoV, % 

3.3 

7.4 

4.4 

4.2 

8.4 

4.4 


TABLE 3. RESULTS FOR THREE-POINT BENDING TESTS. 


Specimen 

Damage 

initiation 

load 

kN 

Damage 
initiation 
flange strain 
pe 

Damage 
initiation 
skin strain 
pe 

Flange 

debond 

load 

kN 

Flange strain 
at flange 
debond load 
pe 

Skin strain 
at flange 
debond load 
pe 

11 

404 

3207 

3811 



MEM 

13 

433 

3051 

3691 

484 


iB 

15 

445 

3231 

3659 

468 


3868 

17 

425 

3036 

3657 

434 



19 

431 

3023 

3481 

488 

3428 

3945 

Mean 

428 

3110 


463 

3370 

3957 

Std. Dev. 

14.9 

101 

118 

24.4 

155 

186 

CoV, % 

3.5 

3.2 

3.2 

5.3 

4.6 

4.7 













TABLE 4. RESULTS FOR ATB TESTS. 


Specimen 

Flange 
debond 
load, kN 

Stroke at 
flange debond 
load, mm 

Flange strain at 
flange debond 
load, pe 

Skin strain at 
flange debond 
load, pe 

12 

2.8 

31.6 

1318 

7199 

14 

- 

- 

- 

- 

16 

2.9 

33.2 

1232 

7254 

18 

2.9 

33.9 

1276 

7295 

20 

2.2 

25.1 

1278 

7015 

Mean 

2.7 

31.0 

1276 

7191 

Std. Dev. 

0.3 

4.0 

35 

124 

CoV, % 

11.6 

13.0 

2.8 

1.7 


axial load = 17.8 kN 


TABLE 5. RESULTS FOR FINITE ELEMENT STRESS ANALYSIS. 


Specimen 

Element 1584 

Mp a 

Element 1604 

MPa 

Element 1624 
°max> MPa 

Element 1 644 

°max> MPa 

Tension 

51.4 

70.2 

87.7 

100.9 

Bending 

46.7 

62.9 

78.0 

89.8 

ATB 

65.0 

93.1 

119.7 

139.0 


for comparison, transverse tensile strength: 61.1 MPa for AS4/3501-6 [12] 


TABLE 6. RESULTS FOR DELAMINATION GROWTH IN 90745° INTERFACE. 


Specimen 

G p J/m 2 

G n , J/m 2 

G x , J/m 2 

gjg t 

G c , J/m 2 

Tension* 

53.4 

13.3 

66.7 

0.200 

118 

Bending* 

42.4 

10.2 

53.0 

0.200 

118 

Tension 1 

63.8 

15.7 

79.5 

0.200 

118 

Bending + 

49.5 

12.4 

61.9 

0.200 

118 

ATB + 

67.3 

24.1 

91.4 

0.260 

132 


* results computed at damage initiation load + results computed at flange debond load 













TABLE 7. RESULTS FOR DELAMINATION GROWTH IN FILM/0 0 INTERFACE. 


Specimen 

Gj, J/m 2 

G n , J/m 2 

G t , J/m 2 

g,/g t 

G c , J/m 2 

Tension* 

362 

76.8 

439 

0.175 

112 

Bending* 

272 

39.0 

311 

0.125 

101 

TensioiT 

467 

93.1 

530 

0.189 

114 

Bending + 

318 

45.9 

364 

0.130 

103 

ATB" 

358 

191 

549 

0.349 

155 


* results computed at damage initiation load + results computed at flange debond load 











Figure 3. Three-Point Bending Test Set-Up. 






I 










Axial Actuator 


Transverse Load Cell 


Transverse Actuator 


Connector Pin 


LVDT 


Specimen 


Bottom Grip 


Figure 5. ATB Test Set-Up Close View. 











25 


Load, kN 


Strain 




Figure 7. Typical Strain-Load Plot for Tension Tests. 
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Figure 8. Typical Load-Displacement Plot for Three-Point Bending Tests. 
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Figure 9. Typical Strain-Load Plot for Three -Point Bending Tests. 
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Figure 10. Typical Load-Displacement Plot for ATB Tests. 
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Figure 11. Typical Strain-Load Plot for ATB Tests. 
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Figure 12. Specimen with Crack Locations. 
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Figure 13. Typical Damage Patterns. 
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(a) Delamination A in the 90745° Flange Ply Interface at Comer 4 
Matrix Cracks Delamination B1 Delamination B2 



(b) Delamination B1 in the Top 0° Skin Ply and Beginning of Delamination B2 
in The top 0°/45° Skin Ply Interface at Comer 2 


Figure 14. Micrograph Showing the Side View of a Failed Three-Point Bending Specimen. 
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Figure 17. Loads and Boundary Conditions For Tension and 
Three-Point Bending Case. 







Figure 18. Loads and Boundary Conditions for ATB Test. 



Figure 19. Virtual Crack Closure Technique (VCCT). 
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(a) Tension Specimen 




(b) Bending Specimen 
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(c) ATB Specimen, Scale Different from (a) and (b) 


Figure 20. Deformed Test Specimen Geometries 
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Figure 21. Comparison of Measured Displacements with Computed Results 

for Tension Load Case. 
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Figure 22. Comparison of Measured Strains with Computed Results 
for Tension Load Case. 
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Figure 23. Comparison of Measured Displacements with Computed Results 
for Three-Point Bending Load Case. 
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Figure 24. Comparison of Measured Strains with Computed Results 
for Three-Point Bending Load Case. 
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Figure 25. Comparison of Measured Displacements with Computed Results 

for ATB Load Case. 
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Figure 26. Comparison of Measured Strains with Computed Results 

for ATB Load Case. 
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Figure 27. Calculation of Force and Moment Resultants. 
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Figure 28. Computed Force and Moment Resultants at Flange Tip. 





Figure 29. Trajectories of Maximum Principal Stresses. 
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Figure 30. Comparison of Principal Tensile Stresses 
Computed at Damage Initiation Load. 
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Figure 31. Computed Strain Energy Release Rates for Delamination Growth 
in a 90745° Flange Ply Interface for Tension Tests. 
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Figure 32. Computed Strain Energy Release Rates for Delamination Growth 
in a 90745° Flange Ply Interface for Three -Point Bending Tests. 
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Figure 33. Computed Strain Energy Release Rates for Delamination Growth 
in a 90745° Flange Ply Interface for ATB Tests. 
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Figure 34. Computed Mixed Mode Ratios for Delamination Growth 
in a 90745° Flange Ply Interface. 
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Figure 35. Computed Strain Energy Release Rates for Delamination 
Growth Between Adhesive and 0° Skin Ply for Tension Tests. 



Figure 36. Computed Strain Energy Release Rates for Delamination 
Growth Between Adhesive and 0° Skin Ply for Three-Point Bending Tests. 
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Figure 39. Mixed-Mode Delamination Criterion for AS4/3501-6 [13]. 
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Figure 40. Comparison of Computed Total Strain Energy Release Rates in a 
90°/45° Flange Ply Interface with Mixed-Mode Fracture Toughnesses. 
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Figure 41. Comparison of Computed Total Strain Energy Release Rates in a 0° Skin 
Ply/ Adhesive Film Interface with Mixed-Mode Fracture Toughnesses. 





Figure 42. Comparison of Computed Total Strain Energy Release Rates in a 90°/45 
Flange Ply Interface with Mixed-ModeFracture Toughnesses. 





Figure 43. Finite Element Model of a Damaged Specimen. 
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Figure 44. Computed Strain Energy Release Rates for Delamination Growth in a 90745° Flange 
Ply Interface with and without a Preexisting Disbond Between Adhesive and 0° Skin Ply. 
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Figure 45. Comparison of Computed Total Strain Energy Release Rates in a 
90745° Flange Ply Interface with Mixed-Mode Fracture Toughnesses. 
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Figure 46. Comparison of Measured Displacements with Computed Results 
for Three-Point Bending Load Case. 
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Figure 47. Comparison of Measured Strains with Computed Results 
for Three-Point Bending Load Case. 
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Figure 48. Comparison of Measured Displacements with Computed Results 

for ATB Load Case. 


8000 
7000 
6000 
5000 
Strain 4000 
3000 
2000 
1000 
0 

0 0.5 1 1.5 2 2.5 3 3.5 4 

Transverse Load, kN 

Figure 49. Comparison of Measured Strains with Computed Results 

for ATB Load Case. 






Strain 
Energy 
Release 
Rate G, 

J/m 2 



Figure 50. Computed Strain Energy Release Rates for Delamination Growth 
in a 90745° Flange Ply Interface for Three-Point Bending Tests. 
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Figure 5 1 . Computed Strain Energy Release Rates for Delamination 
Growth Between Adhesive and 0° Skin Ply for Three-Point Bending Tests. 
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Figure 52. Computed Strain Energy Release Rates for Delamination Growth 
in a 90745° Flange Ply Interface for ATB Tests. 
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Figure 53. Computed Strain Energy Release Rates for Delamination 
Growth Between Adhesive and 0° Skin Ply for ATB Tests. 




Figure 54. Comparison of Computed Total Strain Energy Release Rates in a 
90745° Flange Ply Interface with Mixed-Mode Fracture Toughnesses. 



Figure 55. Comparison of Computed Total Strain Energy Release Rates in a 0° Skin 
Ply/ Adhesive Film Interface with Mixed-Mode Fracture Toughnesses. 
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